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Abstract 

We investigate the twist four longitudinal structure function of deep 

inelastic scattering and show that the integral of — is related to the ex- 
pectation value of the fermionic part of the light-front Hamiltonian density 
at fixed momentum transfer. We show that the new relation, in addition 
to providing physical intuition on relates the quadratic divergences of 

F£^^ to the quark mass correction in light-front QCD and hence provides 
a new pathway for the renormalization of the corresponding twist four op- 
erator. The mixing of quark and gluon operators in QCD naturally leads 
to a twist four longitudinal gluon structure function and to a new sum rule 
J dx^ = 4-^, which is the first sum rule obtained for a twist four observ- 
able. The validity of the sum rule in a non-perturbative context is explicitly 
verified in two-dimensional QCD. 

PACS: ll.lO.Ef, ll.lO.Gh, ll.55.Hx, 13.40.-f 
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1. Introduction 



Measurements [0] of R = ^ show [Q] that dynamical higher twist effects play an im- 
portant role in nucleon structure experiments in the SLAG kinematic range. Indeed, twist 
four ("power suppressed") contributions to the deep inelastic structure functions contain 
non-trivial, non-perturbative information on the structure of hadrons 0. So far, however, 
a clear and simple physical picture of twist four contributions is lacking. In addition, power 
counting indicates that in the bare theory, twist four-matrix elements will be afflicted with 
quadratic divergences 0]. Understanding the nature of these divergences and how to remove 
them (renormalization) is essential in order to compare with the experimental data. 

At twist two level, light-front analysis of deep inelastic scattering provides an intuitive 
physical picture of various structure functions. Recently, the extraction of the unpolarized 
twist two structure function F2 from the transverse component of the bilocal matrix element 
and the physical picture of the transverse polarized structure function gx and associated 
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issues have been analyzed in a light-front field theory approach [|]. The method uses 
Fock space expansion for the Hamiltonian which is physically intuitive and is straightforward 
to calculate. 

To resolve the outstanding issues mentioned above, we use this same framework to eval- 
uate the twist four contribution to the longitudinal structure function. The resulting rela- 
tionship with the fermionic part of the light-front Hamiltonian density provides a convenient 
pathway to renormalization and leads to a definition of a twist four longitudinal gluon struc- 
ture function. We then find a sum rule, the first one at the twist four level. 

The BJL limit, together with light-front current algebra |^ in A'^ = gauge, the tools 
used in the pre-QCD era, lead to the twist four part of the fermionic contribution to the 
longitudinal structure function 
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where / represents a quark (g) or anti-quark (q) or both depending on the target | P), and 
where the twist two contribution to the -p2(/) structure function 
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The bilocal current operator. 
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Note that in the case of quark (anti-quark) contributions, the second (first) term in the 
expression for the bilocal current in Eq. (P) vanishes. We have. 
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Using ijj{y )7 -0(0) = 2ip {y )ip (0), where 

^'{z) = jj dy-e{z- - y') [«^.(z<9^ + gA^) + 7°m] ^'+(2/-), 
P f ^Q^^^ yffQ arrive at 
dy~e--^^^y~''{P I ^p^^\y~) 
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and 



Mo = -2- 
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Thus we have obtained an expression for the twist four part of the fermionic component of 
the longitudinal structure function. After the establishment of QCD as the underlying theory 
of strong interactions, the twist four part of the quark contributions to the longitudinal 
structure function has been given in the limits of vanishing target transverse momentum 



and massless quark and in the = gauge using the QCD factorization method |[10|JTT 
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J dy-e-i^^y-^{P I i^^\y-) 
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Using Eqs. @, and (^), and taking the limit of vanishing target transverse momen- 
tum and massless quark, our result given in Eq. (|l]) reduces to that obtained via the QCD 
factorization method (Eq. p!ol) . 

Note that the expression for F^"^ given in Eq. (|I]) appears to violate transverse boost 
invariance. But, we exhibit below with explicit calculations in Sec. 3 that the P-*- dependence 
cancels between Eqs. (||) and (P) so that the full Flj^^ is indeed boost invariant. 

2. Sum rule for the twist four longitudinal structure function 

From Eq. (H) it follows that F2[f){—x) = F2[f){x) and from Eq. (0) we explicitly find 
that P^(^j(— x) = — P£y-j(a;). Consider the integral 
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Interchanging the orders of x and y integrations and carrying out the integrations explicitly, 
we arrive at []T2[ 



dx- 
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:(P|2V^7+a+^/;|(0)|P)l. (12) 



Identifying iip^'d'^i/j = 6^ , the fermionic part of the light-front QCD Hamiltonian density 
and iip'j'^ ijj = Qq^-, the fermionic part of the light-front QCD longitudinal momentum 
density, (see Eqs. (|15D and ([T7| ) below), we have: 
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where we have used the fact that the physical structure function vanishes for a; > 1. 

The integral of — ^ is therefore related to the hadron matrix element of the (gauge 
invariant) fermionic part of the light-front Hamiltonian density. This result manifests the 
physical content and the non-perturbative nature of the twist-four part of the longitudinal 
structure function. 

The fermionic operator matrix elements appearing in Eq. (|13[) change with as a result 
of the mixing of quark and gluon operators in QCD under renormalization. Next we analyze 
the operator mixing and derive a new sum rule at the twist four level. 

The symmetric, gauge-invariant energy-momentum tensor in QCD is given by 

-g'^'ii (i^^Dx - m) i/j. (14) 

The last term vanishes using the equation of motion. 

Formally, we split the energy momentum tensor into a "fermionic" part 9^'^ representing 
the first term in Eq. ( [T^ ) and a "gauge bosonic" part 9^'^ representing the second and 
third terms in Eq. (|l^. To be consistent with the study of deep inelastic structure function 
formulated in the A'^ = gauge, we shall work in the same gauge. We have, for the fermionic 
part of the longitudinal momentum density, 

9++ = hp^+d+ifj. (15) 

For the fermionic part of the Hamiltonian density, we have, 

9+- = ^^p+^^-1JJ+ + gip+^A-ip^ + i^j-^d^i)-. (16) 

Using the Dirac equation for the fermion, we find that the sum of the first two terms equals 
the third term in the above equation. Therefore, 

9+-= i^^-d^ij = 2iij-^d^ip- (17) 
= 2ij+^\a^.{id^ + c^A^) +7%]^[a^.(za^ + gA^) + -f^mU+. (18) 

The gauge boson part of the Hamiltonian density is given by 



= {d'Aif + 2gr'"^A'^Aid'Ai + A^.A^^A^A^^ 



+ 2gd'A\ {r'^J^,d+A^^ + 2(^+)tT>+) 

+ g^ (^) {r'^^d^K + 2(^+)tr>+) {r'^'A^d+A^, + 2(v^+)trV 



(19) 



where we have used the equation of constraint for the gauge field. 
We define the twist four longitudinal gluon structure function 
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:(P|r+(0) IP)]. (21) 



(P I ^+-(0) I P) = 2P+P- = 2(M2 + (P^)2) and (P | ^++(0) | P) = 2(P+)2, (22) 

where M is the invariant mass of the hadron. Thus we arrive at the new sum rule for the 
twist four part of the longitudinal structure function 



^dx 



(23) 



To our knowledge, this is the first sum rule at the twist four level of deep inelastic scattering 
or for QCD in general. 
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3. Implications of the sum rule 

Next, we investigate the implications of Eq. ( ]T3| ) for quadratic divergences in ^l^y For 
simplicity, we select a dressed quark target and evaluate the structure functions to order g'^. 
That is, we take the state | P) to be a dressed quark consisting of bare states of a quark 
and a quark plus a gluon: 

\P,a)=M\P,^)\0) 



+ E / 7^ \ / ; ^ \ V2(27r)3p+^^(P -k,- k,) 



02(P,a I ki,ai-k2A2)h\ki,ai)a\k2A2) \ 0). (24) 

Without loss of generality, we work in the massless limits for the dressed and bare quarks 
to obtain 



-Mi= ^ / dy-e-^^P^y-%P I r\y-) 



[id^ + gA^{y)].a^a'-.[id^ + gA\Q)]^+{Q) \ P), (25) 



and 



A^2 = -4^xF2(,)(x). (26) 

Note that the matrix elements appearing in Eqs. (§), (^, (0), and ([T9|) involve prod- 
ucts of operators. In this initial study of quadratic divergences we treat them as normal 
ordered which is sufficient for our purposes here. This aspect of the problem merits further 
investigation which we hope to undertake in the near future. However, the terms we drop 
are quadratically divergent and they will affect only the counterterm structure. 

First we evaluate the contribution A^2 given in Eq. (|26|). We obtain, 

iP' 



2 ^, 1 I ^2 ^ 1 I „,2 



-^2 = -ACf^—rJ-x' \6{1 - x) + ^lnA'[^r±^ - 5{l - x) f dy^-±^]\ (27) 
Q^L In 1 — X J 1 — y ^ 

where we have cutoff the transverse momentum integral at A. Note that the result in Eq. 
(p7|) violates transverse boost invariance. 
We have, from Eq. (pSj), 



Mi= / dy-e-'^^'^y-^P I ,p^\y-)id^ fi;+iO) \ P) 

+ I dy~e-^^^''y'%P \ i)^\y-)id^ .a^a^ .A^{^)i)^{Q) \ P) 

+ / dy-e-'^^'^y-^P I i;^\y-)A\y).A\0) \ P) (28) 

= Ml + Ml + Ml + Mf. (29) 
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Since the operators in Eq. (pSj) are taken to be normal ordered, the contribution of A^f 
vanishes to order g^. 

Exphcit calculation leads to the diagonal Fock basis contributions 



{M.Uag = M1 = 4Cf^—fx'\6{l -x) + ^lnA'[-±^ - 5{1 - x) ( dy- 

(J^ L ZTX I — X J 



y 



+ Cfh7V.^- (30) 



g"^ AM + x^ 

The first term here explicitly cancels the term A^2 given in Eq. (|27|) . 
Off-diagonal contributions 

iM,)nonMag = M', + = -Cf^^-^. (31) 



Adding all the contributions, we have, 

'^^''^'--Cf:^.^-^M + -). (32) 



9^ A2 

L(,)l^-J--W^g-2 



As anticipated from power counting, we have generated quadratic divergences for Fl in 
the bare theory. Now we are faced with two issues: (a) What is the principle for adding 
counterterms? (b) What determines the finite parts of the counterterms? 

Further, since Fl^^^ is directly related to the physical longitudinal cross section, we 
expect Fl(q) to be positive definite (see, for example, Ref. |]10|)- From our answers we see 



that we get a negative answer which is free from end point singularities but is quadratically 
divergent. Note that it is easy to show from our expressions in Eqs. and (^ for a free 
quark of mass m, FJJ^ = 4^5(1 — x), a well-known result. 

Now we show that the sum rule relating the integral of — to the fermionic part of 
the light-front Hamiltonian density helps us to understand the results we obtained for 
for a dressed quark in the bare theory and indicates how to add counterterms to renormalize 
the theory. 

A straightforward evaluation leads to 



{P I (0) I P)non<i^a, = 'C f\^-,k' f (33) 

^ In'^ JO x \ — X 



{P I eriO) I F),„, - I ^r(O) I PUa, = Cfl^.A'l' (34) 

Adding the diagonal and off-diagonal contributions from the fermionic part of the Hamil- 
tonian density, we arrive at 

Q,[iP\n-('»\P)-y^iP\^ri'>)lP)]^-C,l,^J^ -(1+.). (36) 

Comparison of Eqs. (^) and (|35|) immediately shows that the relation given in Eq. ([I3D 
has been successfully tested to order g"^ in perturbative QCD in the bare theory. 
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Apart from the verification of the sum rule, the dressed quark calculation serves to 
illustrate the physical origin of quadratic divergences in the corresponding twist four matrix 
element. The sum rule explicitly relates the quadratic divergences in F^^^^ to the fermion 
mass shift in light-front Hamiltonian perturbation theory. The quadratic divergence is to be 
removed by adding counterterms. The precise selection of counterterms is dictated by the 
regularization and renormalization of the light-front QCD Hamiltonian, the diagonalization 
of which results in the hadronic structure. The choice of counterterms in the Hamiltonian, 
in turn, dictates the counterterms to be added to the longitudinal structure function which 
results in the prediction of theoretical twist four contribution to Fl. 

To complete the discussion of the fermion mass shift, we now consider the contributions 
from the gluonic part of the Hamiltonian density. 

Take the off-diagonal contribution arising from 

1 ,+t™ 



M3 = -{P\ 4c/— (aM*'^jV^+'T>+ I P). (36) 



A straightforward evaluation leads to 
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Thus from Eqs. (^3]) and (0) we have, 

{P I K-(0) + (0)] I P)non..a, = 'Cf^,^' ^ (T^' ^^^^ 

Next let us consider diagonal contributions. From the gauge boson part of the Hamilto- 
nian density, we have. 

Adding the diagonal contributions from the fermion and gauge boson parts, i.e., Eqs. (|34D 
and (|39|) we arrive at 



{p I ^^-(0) I pu., - |^(p I r+(o) I pu., = Cf£,A^l' ^(1^- 

Thus the total contribution (quark and gluons) to the expectation value of the Hamiltonian 
density is given by 

( P^V dr 1-1- 

I "^-(0) I P) - y^^P I ^-(0) I P) - -C^AjA^/^ (41) 

This result is directly related to the mass shift of the fermion to order g^ in light-front 
perturbation theory. We have in the limit of zero bare mass, (see Eq. (4.10) in Ref. 
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Note that in the massless hmit we encountered only quadratic divergences in the twist 
four part of the longitudinal structure function. It is now obvious that for a massive quark, 
we will also encounter logarithmic divergences. 

To test the sum rule given in Eq. (|23| ) explicitly in a non-perturbative context we turn 
to two-dimensional QCD. In 1+1 dimensions, in A'^ = gauge, we have, 

1 r/r 9 r n 

-Fl^ix) = 7^{P I e-(0) + (0) I P), (43) 



x-^(«)^-^^-Q^ 

with e+- = 2mV^^ -Jt^"*" and 9+- = -4:g^ij+^T''tjj+jQ^tjj+^T''^+. In the standard one 
pair (qq) approximation to the ground state, with the help of the bound state equation ('t 
Hooft equation) obeyed by the ground state wave function i/j^x) for the meson 

MV(.r) = - / dy *^fl^ (44) 

x{l — X) ' IT J [x — yY 

together with the normalization condition /q^ dxil)'^ (a;) = 1 , it is easily verified that the twist 
four longitudinal structure function of the meson obeys the sum rule 

^f^fr.-l(P|r-(0)|P>^4^. (45) 

To summarize, we have investigated the twist four contributions to the longitudinal 
structure function in deep inelastic scattering in an approach based on Fock space expansion 
methods in light-front field theory. Our motivation has been to gain physical intuition on 
the twist four part of the longitudinal structure function and to understand the occurrence 
of quadratic divergences in FJJ^^ and the associated renormahzation issues. We have found 

that the integral of — is related to the expectation value of the fermionic part of the 
light-front Hamiltonian density. The relation, which is tested to order g'^ in perturbative 
QCD in the bare theory, explicitly relates the quadratic divergences generated in F^^^ to 
the fermion mass shift arising from the fermionic part of the light-front Hamiltonian density. 
By investigating the operator mixing problem, we have derived a new sum rule which relates 

the integral of ^ ^ to the invariant mass of the hadron. We have explicitly checked the 
validity of this sum rule in two dimensional QCD. 

The phenomenological consequences of the new sum rule are worth studying. Clearly, 
experiments to measure the twist four longitudinal structure function will reveal the fraction 
of the hadron mass carried by the charged parton components. These experiments will 
complement the longitudinal momentum and helicity distribution information obtained at 
the twist two level. It is of interest to investigate the feasibility of the direct measurement 
of the twist four gluon structure function in high energy experiments. Recent work of 
Qiu, Sterman and collaborators have shown that semi-inclusive single jet production in 
deep inelastic scattering Ijl^ and direct photon production in hadron nucleus scattering |l^ 



provide direct measurement of twist four gluon matrix elements. 

On the theoretical side we note that some significant progress has been made recently in 
the bound state problem in light-front QCD [|TB[. In the near future, we plan to undertake a 



non-perturbative calculation (utilizing Fock space expansion and Hamiltonian renormahza- 
tion techniques) of the longitudinal structure function for a meson-like bound state. Such a 
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calculation will undoubtedly help to check the validity of current phenomenological models 
based on simple assumptions Jl^ employed in analyzing the data. Whether the relation 

between the integral of — and the hadron expectation value of the fermionic part of 
the light-front QCD Hamiltonian density can provide insight in to the small x behavior of 
^I'iq) is ^ls° worth investigating. 

One of the authors (AH) would like to thank Wei Min Zhang for useful discussions. This 
work is supported in part by the U.S. Department of Energy under Grant No. DEFG02 - 
87ER40371, Division of High Energy and Nuclear Physics. 
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